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Abstract 

A unified theory of material defects, incorporating both the smooth and the 
singular descriptions, is presented based upon the theory of currents of Georges de 
Rham. The fundamental geometric entity of discourse is assumed to be represented 
by a single differential form or current, whose boundary is identified with the defect 
itself. The possibility of defining a less restrictive dislocation structure is explored in 
terms of a plausible weak formulation of the theorem of Frobenius. Several examples 
are presented and discussed. 
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Introduction 

While the precise definition of the concept of material defect must be left to every par- 
ticular context, a common feature of all theories dealing with defects (dislocations, inho- 
mogeneity, and so on) appears to be that the presence of defects translates itself mathe- 
matically into the lack of integrability of some geometric entity. In a general differential 
geometric framework, questions of integrability pertain to differential forms and their ex- 
actness or lack thereof. It seems appropriate, therefore, to undertake a unified treatment 
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of defects by associating to any possible structure under consideration one or more dif- 
ferential forms. On the other hand, since differential forms are, by definition, smooth 
entities, it would appear that the rich variety of isolated defects, whose practical and 
historical importance cannot be denied, might be left out and that a unified treatment 
encompassing both the continuous and discrete cases would remain out of reach of a sin- 
gle formal apparatus. The situation is similar in many other engineering applications, 
where concentrated entities (forces, masses, charges) can be seen as limiting cases of their 
smooth counterparts. The unified mathematical treatment of these cases was historically 
achieved by the theory of distributions, where the singular entities are represented not 
by functions but rather by linear functionals on a suitable space of test functions. The 
most common example is provided by the Dirac delta which assigns to each compactly 
supported function in R its value at the origin. Since a scalar field can be considered as 
a particular case of a differential form, it is not surprising that L. Schwartz's theory of 
distributions can be extended to forms of all orders. This extension, achieved by G. de 
Rham [?], is completely general and independent of any metric considerations, a feature 
that should be considered essential in a truly general geometric setting. De Rham in- 
troduced the terminology of currents to designate his generalized differential forms. It is 
this tool that will serve our purposes in the present formulation of the unified theory of 
defects!*] 

1 Currents 
1.1 Definition 

A p-current on an n-dimensional manifold is a continuouJl] linear functional T[(j)\ on 
the vector space of all forms with compact support in A^. To understand in what 

sense this definition is consistent with that of smooth forms, it suffices to exhibit the 
latter as a particular case of the former. Let, therefore, u represent a smooth p-form on 
A4. We can uniquely associate to it the (n — p)-current defined as the linear operator 



for all (n — p)-forms (p with compact support in Ai. Strictly speaking, the [n — p)- 
current cannot be "equal" to the p-form u, but they are indistinguishable from each 

*This work extends our previous article [?]. 

'I'By continuity we mean that the sequence of evaluations T[(f)i\ on a sequence of C°° p-forms supported 
within a common compact subset of Ai tends to zero whenever the coefficients and all their derivatives 
of a coordinate representation of the forms (f>i tend to zero uniformly as i — >■ oo. 




(1.1) 
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other in terms of their integral effect on all "test forms" 0. Thus, a form bears to its 
associated current the same relation that a function bears to its associated distribution. 
An important non-trivial example of a current that is not associated to any differential 
form is the following. Let s be a p-simplex in A4. We associate to it the p-current defined 
by: 

Ts[(l>] = !<!>, (1.2) 

s 

for all compactly supported p-forms (p. The definition above can be extended by linearity 
to arbitrary chains. These examples show how an integrand and a domain of integration 
are unified under the single formal umbrella of currents. 



1.2 Operations 

(1) Currents of the same dimension can be added together and multiplied by real num- 
bers in an obvious way. 

(2) The product of a p-current T with a g-form a is the (p — g)-current T\_a defined as: 

{TLa)[(j)] = T[a A (j)]. (1.3) 

Similarly, 

ajT = (-l)("-f)9 Tlck. (1.4) 

(3) The product of a p-current T with a vector field X is the (p -|- l)-current: 

{TAXM = T{X^4>]- (1.5) 

(4) The boundary of a p-current T is the (p — l)-current 

ar[0] = T[c/0]. (1.6) 

Using Stokes' theorem for chains, it is easy to show that, for any chain c, 

or, = Tae. (1.7) 



1.3 Possibilities 

The notion of currents opens the doors for generalizing classically smooth differential 
geometric objects, such as connection, torsion and curvature. While this idea is beyond 
the scope of this paper, it is not difficult to intuit the possibilities. Consider, for example, a 
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(non-zero) decomposable differential ]3-form a; on an n dimensional differentiable manifold 
Ai. Thus, there exist p linearly independent 1-forms Ui [i = 1, ...,p) such that: 



uj = Ui A ... A Up. 



(1.8) 



Such a form uniquely determines at each point of x G an (ri — p)- dimensional subspace 
Hx of the tangent space T^Ai. A vector v G T^M belongs to if uji{v) = for each i = 
1, ...,p. The collection Ti of all the subspaces is called a (geometric) (n—p)- dimensional 
distribution on Ai . Conversely, given a distribution, the corresponding decomposable form 
u is determined up to multiplication by a scalar field a : — t- M. A submanifold S of A4 
is called an integral manifold of the distribution "H if for every s G 5 we have TsS = Hg. 
A distribution is completely integrable if at every point x it admits an integral manifold 
of maximal dimension (i.e., n — p). According to (one of the versions of) the theorem of 
Frobenius, a distribution 7i defined by a decomposable form u is completely integrable 
if, and only if, there exists a 1-form /3 such that: 



This is tantamount to saying that, for some choice of the scalar degree of freedom a, 
the form au is closed, namely, there exists an integrating factor a such that d{au) = 0. 
So far, we have been dealing with the smooth case. Assume now that we have a means 
of characterizing the decomposability of a current (perhaps as the limit of a sequence of 
non-zero decomposable forms). We could now declare that a decomposable p-current T 
determines a p-dimensional singular geometric distribution on Ai and define the complete 
integrability of the singular distribution by the condition: 



for some 1-form A stronger condition would be to require that /3 be closed. Notice 
that since connections in general can be regarded as (horizontal) distributions on fibre 
bundles, and since the curvature of a connection is related to its complete integrability, we 
can expect that singular connections can be introduced by means of decomposable currents 
and their non-vanishing curvature can be detected by the violation of a condition such as 
fll.lOp . By this means, a situation is envisioned in which the standard curvature vanishes 
almost everywhere and is concentrated, as it were, at a single point. We remark that 
Equation (11. 101) is by no means the result of a theorem, but only a possible definition of 
complete integrability of a singular distribution. Clearly, the fact that /9 is a smooth form 
may severely limit the singular distributions that can be considered completely integrable. 

"''As a curiosity, it is interesting to remark that Equation (jl.lOl) can be informally regarded as the 
eigenvalue problem of the boundary operator d. Its "eigenvectors" are the completely integrable currents. 



doj = (3 Au. 




dT = (3^T, 
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2 Bravais hyperplanes 



2.1 The smooth case 

The traditional heuristic argument to introduce continuous distributions of dislocations 
in crystalline materials calls for the specification of a frame field (or repere mobile) in the 



obtained by means of a co-frame field, which can be regarded as an M"'-valued 1-form 
on Ai. This point of view suggests perhaps that n linearly independent 1-forms might 
constitute a convenient point of departure for our desired generalization. Each of these 
1-forms would represent a family of Bravais planes. It comes as a surprise, however, that 
defects are meaningful and detectable with just a single family of such planes or, more 
specifically and less surprisingly, that integrability conditions can be associated with a 
single 1-form on a manifold. Geometrically, a 1-form (always decomposable) induces an 
{n — l)-dimensional distribution, that is, a field of hyperplanes. It is physically important 
to point out that, relinquishing the multiplicative degree of freedom alluded to in the 
previous section, a 1-form also specifies a local density of these Bravais hyperplanes. 
Indeed, a covector acting on a vector space defines a family of parallel hyperplanes and 
the evaluation of the covector on a given vector can be pictorially regarded as the 'number 
of hyperplanes' pierced by the vector. We have at our disposal, therefore, two somewhat 
different images induced by the specification of a 1-form a; on a manifold Ai. The first one 
would look at the integrability of the form itself by demanding that u be closed, namely: 



Physically, this condition means not only that the induced distribution is completely 
integrable, but also that the local hyperplane densities are mutually compatible. More to 
the point, we expect the Bravais hyperplanes to fit well with their neighbours not only as 
hyperplanes but also as stacks thereof at each point. The second image, on the other hand, 
would demand only the nice fit between the hyperplanes themselves. The corresponding 
(less demanding) integrability condition would take the form of a Frobenius condition. In 
either case we will refer to the form a; as a layering form. 

§This feature is present also, albeit with the degree of freedom afforded by material symmetries, in 
the constitutively based approach propounded by Kondo [?] and Noll [?] , whereby points are compared, 
in a groupoid-like fashion, via material isomorphisms between their tangent spaces. 



body manifold, and the consequent distant parallelism!^ An alternative, dual, pict 




ure IS 



du = 0. 
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2.2 Dislocations and currents 



From the smooth case we learn two facts: (1) A basic entity to be analyzed for defective- 
ness is expressed in terms of a differential form u. In the case of a 1-form, the entity is a 
distribution of Bravais hyperplanes with their corresponding stacking densities. (2) The 
defectiveness of the structure is measured by the exterior derivative Vt = du oi the form 
representing the basic entity. It is natural, therefore, to identify VL with the dislocation 
or, in more general terms, with the imperfection. Notice that different entities u may 
have the same boundary (if they differ by a closed form), which means that the same 
dislocation structure may arise from different physical objects. 

The generalization of these notions to the non-smooth case is straightforward, provided 
one bears in mind de Rham's invention. Let a p-current T represent some basic physical 
object. Then we call its boundary D = dT the associated dislocation current. Notice 
that D is a (p — l)-current. We say that the object represented by T is defect-free if the 
integrability condition 

dT = 0, (2.2) 

is satisfied. In other words, T is defect-free if it has a vanishing boundary D = 0. An 
important non-trivial example is provided by a p-simplex s embedded in the manifold A4. 
As we have seen, we can associate to s the current defined in Equation (11. 2p . For the 
important particular case p = n — 1, the current Tg can be regarded as the specification 
of a Bravais structure concentrated on s rather than distributed over the whole body. 
Physically, the simplex s may be regarded as a cut inside the body where a putative layer 
of atoms has been inserted or removed. Let us calculate the corresponding dislocation 
current Dg. For every compactly supported {p — l)-form 0, we have: 



DM = dTM = Ts[d4>] = 




s 



Thus, we obtain the important result that the dislocation coincides with the current 
associated with the boundary of the embedded simplex. 

2.3 An edge dislocation 

Consider the open cube Ai = (—1,1)^ in with coordinates x,y,z. Let h denote its 
intersection with the (oriented) lower half-plane x = 0, z < 0. We associate with h the 
current: 




h 
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where denotes an arbitrary compactly supported 2-form in . Notice that the supports 
of these forms must be made entirely from (interior) points of the open cube. As a 
consequence of this observation, we obtain: 

- dniij] = Jdilj = Ji; = niiP], (2.5) 

h L 

where L is the open interval (—1, 1) on the y-axis. Thus wc recover the classical textbook 
description of an edge dislocation as the result of the removal of an atomic half plane. 



2.4 An open book 

An interesting example is provided by the 1-form (f) = dO defined in F = ]R^\{0}, that 
is, the real plane devoid of the origin. We denote by p, 9 the usual polar "coordinates" . 
Clearly, these are not legitimate global coordinates for F since (p, 9) and (p, 9 + 27r) 
represent the same point. Nevertheless, the notation d9 is standard and reinforces the 
fact that the 1-form is closed (namely, dcj) = 0) though not exact. If we should propose 
as a layering form on F, the corresponding distribution would look like the set of all 
rays emanating from (but not including) the origin {0}. The corresponding dislocation 
current is D — dcf) — 0, which means that, as far as the set F is concerned, the given 
layering is defect free. We are interested, however, in extending the form (f) to include the 
missing origin. To this end, we define the following associated current on R^: 

T^[a] = J (pAa. (2.6) 

F 

The subtle point in this definition is that the 1-forms a have compact support in rather 
than in F. We are interested to obtain the corresponding dislocation current D, that is: 

D[f3] = dT[f3] = T[d/3] = J <f>Ad(3 = J d{f3(j)), (2.7) 

F F 

where f3 is any zero-form (function) with compact support in and where dcj) = was 
used. Since we cannot use Stokes' theorem directly, wc resort to evaluate the last integral 
over the domain F^ obtained by subtracting from the closed ball of radius e with centre 
at the origin and then going to the limit as e ^ 0. We obtain: 

j dm = lim J dm = lim J /3<f>^ 27r/3({0}). (2.8) 

F Fe dFe 

Thus, the dislocation current is given by Dirac's delta. To obtain a three-dimensional 
version of the above, we consider the same form (f) — d9 in a cylindrical coordinate system 
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p, 0, z and define the domain F as minus the entire ^-axis. This layering form can 
now be pictured as the pages of an open book evenly spread with the spine occupying the 
^-axis. The corresponding dislocation current is now given by: 

oo 

D[^] = dT^[^] = 27r y 73(0, 0, z)dz, (2.9) 
—00 

where 7 = ^idx + ^2dy + ^^dz is a 1-form with compact support in R^. 
2.5 A screw dislocation 

Two currents that differ by a closed (i.e., zero-boundary) current have the same boundary. 
In the context of defects, we may say that two layering currents that differ by defect-free 
current must exhibit exactly the same defects. This observation can have unexpected 
physical interpretations. Indeed, let us consider the open-book layering current just 
introduced and let us define the 2-current: 

S^T^ + aTa,, (2.10) 

where a is a constant. Since dz is a well-defined closed form over R^, so is the associated 
current T^z- Consequently, the boundary of S coincides with the boundary of T^. The 
layering structure corresponding to S consist of applying to the previous 'pages' a uniform 
twist about the 2;-axis. To ascertain that this is indeed the case, notice that, within the 
domain F, S can be regarded as the 1-form = d9+adz, a closed form. Locally, therefore, 
we can write: 

(l)^d(e + az). (2.11) 
In other words, locally the submanifolds with equation: 

9 + az — constant (2-12) 

are integral submanifolds of the distribution generated by (p. These submanifolds describe 
helicoidal surfaces climbing around the z-axis. This screw layering has a dislocation 
current identical to that of the open book. 

3 Frank's rules 

Within the classical theory of dislocations in crystals, a prominent role is played by the 
Burgers vector concept. Dislocations are assumed to occur along lines only. A Burgers 
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circuit in an atomic lattice consists of a quadrilateral path situated on an atomic 'plane' 
transverse to the dislocation line and with equal numbers of atomic cell steps on opposite 
sides. In a perfect crystal, these paths naturally close. The lack of closure (namely, 
the Burgers vector, denoted by b), on the other hand, is interpreted as the presence of 
a dislocation. When, for example, b is parallel to the dislocation line, we have a pure 
screw dislocation. Clearly, as one advances over the dislocation line, the Burgers vector 
may change in magnitude and direction, so that the question arises as to whether this 
change can be arbitrary. Moreover, dislocation lines may meet and branch out, so that 
a similar question arises in these more involved cases. In an important article [?], F. 
C. Frank introduced the notion of the law of conservation of Burgers vectors, formally 
analogous to Kirchhoff 's laws for electrical circuits (charge conservation) or similar laws 
for fluid flow in pipes (mass conservation). As a consequence of this law, several rules can 
be deduced. For example, the Burgers vector along a dislocation line must be constant. 
Moreover, a dislocation line may not end within the crystal, but only at its boundary. 
At a bifurcation, the vector of the entrant trunk is equal to the sum of the vectors of 
the outgoing branches. Given the importance of these rules in applications, we want to 
place them rigorously within the context of the geometrical theory. It will turn out that 
Frank's rules are various expressions of the general topological criterion that establishes 
that the boundary operator is nilpotent of degree 2, that is, the boundary of a boundary 
necessarily vanishes. 

3.1 The constancy rule 

Since Frank's rules deal always with dislocation lines, we need flrst to establish the notion 
of the support of a current [?]. A current T is equal to zero in an open set U if T[0] = 
for all forms compactly supported in U. The support of T is defined as the complement 
of the maximal open set in which T — 0. Accordingly, we say that a dislocation current 
D is a dislocation line if its support is a curve. Considering, for specificity, a (three- 
dimensional) body A4, we investigate the possibility of existence of a curve L within the 
body, whose ends are not points of A4, with the following properties: (1) L is the support 
of a dislocation 1-current D; (2) D is of the form T^l for some real valued function 
u : M ^ R, namely: 



with some abuse of notation (in the sense that, within the integral, the 1-form a represents 
the restriction to L of the original 1-form a with compact support in A4). We will prove 
that both conditions cannot be satisfied simultaneously unless the scalar function u is 




(3.1) 



L 
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actually constant on L. The proof starts by remarking that, as a dislocation current, D 
must be the boundary of some (layering) current S: 



D^T^L^ dS. (3.2) 
Applying the boundary operator, we obtain 

dD = dT^L = ddS = 0. (3.3) 
Evaluating over a zero-form / yields 

dT^Lif] = TulW] ^ Juf^ J d{uf) - J fdu^- J fdu^ 0. (3.4) 

L L L L 

Since / is arbitrary, we conclude that u — constant on L. Thus, if we interpret u as the 
strength of the dislocation, it follows from this proof that the strength of a line dislocation 
must be constant. Recall that we are analyzing the dislocations associated with a single 
layering system, whence the scalar nature of the line dislocation strength. 



3.2 Branching 

If k lines {Li, i = 1, k > 3) meet at a body point X e A1, and if each line is the support 
of a hne dislocation, we have a case of branching. We assume these hues to originate at 
X and to emerge at the topological boundary of M. (which is to be seen, as standard 
continuum mechanics prescribes, as an ordinary differentiable manifold, not as a manifold 
with boundary). It follows that the boundary of each of the given hues Lj consists of 
the single point X. Denoting by the (constant) strength of the dislocation supported 
by Li, we associate to the total system the dislocation current D — J^^^i'^Li- For an 

i 

arbitrary compactly supported zero-form /, considering that D itself must be a boundary 
of a (layering) 1-form, we obtain: 

= dD[f] ^J2aijdf^ f(X) J2 (3-5) 

i J i 

Since / is arbitrary, we conclude that: 

J^ai^O, (3.6) 

i 

which is Prank's branching rule within the scalar context. 
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4 The second integrability criterion 



We indicated in section [2711 that, given a layering form u, there are two different questions 
that one may try to answer, each one leading to a different integrability criterion. The first 
question, which we have been exclusively addressing so far, is whether or not the layers and 
their respective stacking densities fit well together. The general answer to this question 
is provided by the closedness of the layering form, namely du = 0, or, in the singular 
case, the closedness of the layering current, dT = 0. The second question, which we have 
described only in the case of a decomposablj^ layering form u, is whether or not the 
associated distribution is completely integrable. According to the theorem of Frobenius, 
the pertinent condition is the existence of a 1-form (3 such that du = (3 /\ u. Clearly, 
this criterion of integrability is less demanding than the first. Correspondingly, every 
defect-free decomposable layering form oj according to the first criterion is also defect- 
free according to the second, but the converse is not true. Examples are not difficult to 
construct. In fact, every 1-form in given by the expression u = /(x, y)dx + g{x, y)dy + 
dz, for any given smooth functions / and g, gives rise to a completely integrable two- 
dimensional distribution, although u is closed only when the cross derivatives f^y and g^x 
are identical to each other. 

4.1 Coherence at interfaces 

Working in with natural coordinates x,y,z, let S denote the plane z = 0. Let, 
moreover, the upper {z > 0) and lower {z < 0) half-spaces be denoted by and H~, 
respectively. Consider 1-forms = f^dx + f2dy + f^dz and /~ = f^dx + f2dy + f^dz 
smoothly defined on and , respectively, and define a 2-current T by 



for arbitrary 2-forms (j) compactly supported in M^. We regard this form as defining a 
singular 2-dimensional distribution on M^. Its boundary is the 1-current: 




(4.1) 




(4.2) 



^Note that l-fomis are always decomposable. 
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acting on 1-forms with compact support in R^. In this equation, [[•]] denotes the jump 
operator. Assume now that, for this particular layering current T, we want to establish 
the absence of defects according to our first criterion. Setting 

= (4.3) 

identically, we may first choose arbitrary forms ip whose support does not intersect S and 
obtain the conditions: 

df+ = 0, df- = 0. (4.4) 

What these conditions mean is that, as far as the individual layering forms and /~ are 
concerned, there are no defects at interior points of iJ"*" or H and all possible remaining 
dislocations are concentrated on the surface E. Considering now arbitrary 1-forms if) 
whose support does intersect S, we obtain the following extra point-wise conditions on E: 

[[/i]] = 0, m = 0. (4.5) 

The jump [[/a]] of the 2;-component of the layering form can be prescribed arbitrarily. 
The geometric interpretation of these results can be gathered by first considering the case 
in which /f , and f2 vanish altogether while — A and — B, where A and 
B are different constants. We have then a purely horizontal layering (i.e., parallel to E) 
that undergoes an abrupt change of density across E. There are no defects in this kind of 
layering. On the other hand, if fi, fi, and were to vanish identically while f2 — A 
and f2 — B, the layering would be vertical (i.e., perpendicular to E) and there would be 
an incoherence defect across E, unless A = B. In a more general case, we would have 
that the layering consists of leaves that may have a kink upon crossing the surface E, but 
are otherwise continuous. 

We have considered so far the absence of defects according to the first integrabihty crite- 
rion, namely, by checking that the layering current T is closed. We investigate now the 
consequences of demanding only the existence of a 1-form (3 such that 



dT = pjT. (4.6) 



More explicitly. 



Jm^i^ + J dtA^+ J df-AiP 



S H+ H- 

= j /+A(/5A^) + j /-A(/3Az^). (4.7) 

H+ H- 

Considering first 1-forms whose support does not intersect E, we obtain now the con- 
ditions 

df+ = f+Af5, df+ = f+A/3. (4.8) 
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In other words, the distributions induced by /"*" and /~ in the interior of their respective 
domains are required only to be completely integrable. By considering 1-forms '0 whose 
support intersects S, however, we recover the coherence condition (14. 5p . We observe not 
only that this condition could perhaps be further relaxed, but also that, even in terms 
of the complete integrability of the upper and lower distributions, a single form (5 should 
act as integrating factor for both. Clearly, a still weaker form of the Frobenius condition 
f ll.lOP could be postulated in terms of a current multiplier rather than a smooth form (5. 
This issue is the subject of further study. 



4.2 Broken leaves 

We define the 2-current 



r[0] = j ct>+ j aA(l), (4.9) 

S in)3 



where E is the same as in the previous example and a is a 1-form defined over R^. Its 
boundary is: 

D[ip] = dT[tp] = j dilj + jaAdij = JdaAip, (4.10) 



d3 



for arbitrary 1-forms ip with compact support in M . The stronger integrability criterion 
dT = yields the expected condition 

da = 0. (4.11) 

We note that the integral over E is automatically closed as a current and, consequently, 
has no effect on the result when applying the (linear) boundary operator. On the other 
hand, demanding only the satisfaction of Equation fll.lOp . we conclude that a 1-form (3 
must exist such that 

dT[ij] = (/3jT)[^] = J /3AiIj + j aA^AiJj. (4.12) 

This identity implies: 

da = aA/3 on M^ (4.13) 

and 

(3 = onS. (4.14) 

To grasp the meaning of these integrability conditions, we observe that if we were to 
ignore the integral over S in the definition of T in Equation (14.90 . the satisfaction of the 
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integrability condition fl4.13p would imply that the distribution is completely integrable, 
thus constituting a regular foliation of M^. The presence of the integral over S has 
the effect of breaking the leaves. The vanishing of (3 repairs the damage. An intuitive 
realization of this picture can be obtained by considering the following sequence of forms: 

ai = {I - Ci)dy + Cidz, (4.15) 

where Cj = Ci{z) are scalar functions supported on the interval [— 2~\2~*] and such that 

adz = 1. (4.16) 



The sequence of 2-currents 

r,[0] = I a,A(P (4.17) 



converges to a 2-current T whose leaves are of the broken kind described above. 
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